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Evidence for the quasi one-dimensional (ID) antiferromagnetism of CuO is presented in a frame- 
work of Heisenberg model. We have obtained an experimental absolute value of the paramagnetic 
spin susceptibility of CuO by subtracting the orbital susceptibility separately from the total sus- 
ceptibility through the 63 Cu NMR shift measurement, and compared directly with the theoretical 
predictions. The result is best described by a ID S = 1/2 antiferromagnetic Heisenberg (AFH) 
model, supporting the speculation invoked by earlier authors. We also present a semi-quantitative 
reason why CuO, seemingly of 3D structure, is unexpectedly a quasi ID antiferromagnet. 



PACS numbers: 74.25.Ha, 74.72. Jt, 76.60.Cq 



I. INTRODUCTION 

CuO is one of the materials closely related to the 
high-T c cuprate superconductors, especially in light of 
the strong antiferromagnetic correlation in the Cu-O-Cu 
bonds. It may be surprising that unusual magnetic prop- 
erties are found in CuO which consists only of the essence 
of cuprates, copper and oxygen ions, and has a nearly 
3D structure from the viewpoint of chemical bonding. A 
close inspection of the magnetic properties in CuO may 
give us an insight into the underlying physics in the mag- 
netism and superconductivity of the cuprate family. 

Extensive studies including magnetic susceptibility, 
specific heat, photoemission, neutron scattering and 
NMR measurements on CuO have revealed that (1) suc- 
cessive magnetic transitions are occurred at Tjvi = 212 
K and T N2 = 231 K 0, 0, (2) a Neel state with 
the easy axis of [010] direction H is achieved be- 
low Tjvi = 212 K with the superexchange coupling of 
J = 67±20 meV IE Hi which is an order of magni- 
tude larger than that expected from the Neel temper- 
ature, suggesting a strongly correlated and low dimen- 
sional spin system H, Ej , (3) the significantly reduced 
Cu 2+ spin moment 0.65/ib compared with 1/xb expected 
for a Cu 2 + ion is observed 0, U at T = 4 K due 
either to quantum spin fluctuations in low dimensional 
system or to covalent effect, (4) the temperature depen- 
dence of paramagnetic susceptibility 0, 0, IT^ | 
shows a broad peak at around 540 K [l(| , reminiscent of 
quasi ID or 2D antiferromagnet, and (5) CuO belongs to 
the charge-transfer gap insulators [l5j . 

All the 3d transition metal monoxides, with the only 
exception of CuO, are 3D antiferromagnets. It may be 
unexpected to find a low dimensional magnetism in such 
a chemically 3D structure as the monoxide CuO. In fact, 
many experimental data of CuO have been explained in 
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the scheme of quasi ID antiferromagnet. It has been as- 
sumed that the strongest antiferromagnetic coupling may 
reside on a particular Cu-O-Cu bond 0E3 by invoking 
the Anderson model of superexchange interactions which 
tells the larger bond anglepreferred for the stronger an- 
tiferromagnetic coupling [To , IT^ . 

The proposed picture of spin ID chain, however, would 
not be obvious, because the intra-chain bond is even 
longer than the inter-chain ones. We note a close com- 
petition present among the bond angles as well as bond 
lengths in CuO. The ID picture may be conceivable but 
still only a speculation unless quantitative evidence for 
the bond angle scheme is presented. 

We will find below that making a distinction between 
ID and 2D involves a delicate problem. No investiga- 
tion, to our knowledge, has been performed yet to make 
a comparison between ID and 2D in CuO. The following 
questions are still open: (1) Which is the better model to 
describe the paramagnetic state of CuO, ID or 2D? (2) 
If it turns out to be a ID antiferromagnet as has been 
believed so far, what is the cause for a particular bond 
to be magnetically active and others less active? 

In this paper, we present evidence for a quasi ID anti- 
ferromagnetism in CuO beyond 2D antiferromagnetism. 
We carried out the 63 Cu-NMR shift measurement to ob- 
tain an absolute value of the spin susceptibility which 
provides us an opportunity to make direct comparison 
with theoretical predictions. The theories we refer to in 
the present work are those studied in a framework of the 
Heisenberg model. Although the earlier studies about 
susceptibility have shown an experimental indication be- 
ing consistent with the prediction of ID AFH, a lack of 
the knowledge about absolute value of spin susceptibil- 
ity prevents quantitative comparison between experiment 
and theory. We will find below that the absolute value 
is crucial to make a distinction between the ID and 2D 
antiferromagnet . 

A susceptibility measurement gives us only a total sus- 
ceptibility which consists of the three contributions from 
the spin, orbital and ionic-core. Only the spin contribu- 
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FIG. 1: Crystal structure of CuO. The big and small balls 
are copper and oxygen ions, respectively. The left panel: We 
define Cu(l) through Cu(7) for the convenience of discus- 
sion, although all the Cu sites are crystallographically equiv- 
alent. Similarly, oxygen ions, which occupy the other equiv- 
alent sites, are labeled as O(l) through 0(6). The primitive 
directions of [100], [010] and [001] are along lines connecting 
Cu(3) to Cu(6), Cu(5) to Cu(3), and Cu(3) to Cu(7), respec- 
tively. The thin bonds running vertically are along the chain 
direction [10—1], the fat bonds connect the chains. The right 
panel: There are four bond lengths in the crystal. We call 
them a, (3, 7 and 8 in the order of short. The values of them 
are given in Table U] The intra-chain bonds consist of the al- 
ternative sequence of a and 5, the inter-chain bonds similarly 
of (3 and 7. 



tion Xspin is required to be compared with a quantum 
spin theory. A direct measurement of spin contribution 
can in principle be made by neutron experiment, but 
the measurement would be difficult because of poor sig- 
nal intensity in case of weak paramagnetism like CuO. 
Ionic-core contribution Xcore can be estimated from the 
ionic data-base or theoretical calculation. Orbital contri- 
bution Xorb arises from second-order excitations between 
crystal field levels, being material dependent and thus to 
be measured by an experiment like NMR. 

In the absence of experimental information of orbital 
contribution, as is common, its contribution would be 
treated as an adjustable parameter to make a qualita- 
tive comparison in the temperature dependence of spin 
contribution between experiment and theory. The aim of 
our study is to present the experimental absolute values 
of the orbital and spin contributions in CuO to make a 
quantitative comparison possible. 



II. CRYSTAL STRUCTURE AND SAMPLE 
PREPARATION 



TABLE I: The bond lengths (inter-atomic distances) are 
shown in the units of A. We denote the four shortest bonds 
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TABLE II: The bond length L (measured along the bond Cu- 
O-Cu) and bond angle 9 are shown. There are six kinds of 
Cu-O-Cu bond in the crystal. We define the abbreviation 
Bijfe representing the bond Cu(i)-0(j)-Cu(fc). 
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monoxides, rather than the cubic rock salt structure of 
the most. The unit cell contains four copper ions which 
are all crystallographically equivalent, as well as four 
equivalent oxygens. 

A significant feature can be seen from Table [H] and 
Fig. H that the bond Cu(l)-0(2)-Cu(2) (abbreviated as 
B122 in the following) has a noticeably larger bond angle 
compared with all others. This is the reason why this 
bond has been assumed to carry the strongest superex- 
change interaction 0, Il0| . This assumption seems con- 
ceivable, but still only a speculation which needs evidence 
to specify the bond angle and bond length dependences 
of the superexchange interaction. 

We began our sample preparation by initializing the 
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CuO crystallizes in a monoclinic structure with a space 
group of Cc The crystal structure, shown in Fig. 
and Table |H is a unique one among the transition metal 



FIG. 2: The bond length L is plotted against the bond angle 
9. The vertical bars labeled by Bijk denote the bonds Cu(i)- 
0(j)-Cu(fc) as defined in Table El 
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FIG. 3: Polycrystal susceptibility of CuO . The present sam- 
ple (A) and literature data (B) taken from Ref. Ilfl are shown. 
The arrows indicated by T N1 (212 K) and T N2 (231 K) are 
the fst and 2nd Neel temperature, respectively 

stoichiometry of commercially obtained CuO, in which 
the powder sample of the nominal CuO was annealed at 
540 °C in air for a day. We have found that the anneal- 
ing process is crucial to have the Cu NQR/NMR signal 
visible at any temperatures including 4 K and 300 K, 
provably because a possible oxygen deficiency may be 
removed by the annealing process. 

The annealed sample was pulverized to be fine powder 
with 5 (im in typical diameter, and then to be mixed 
with epoxy resin and finally exposed in a magnetic field 
(7 T) to prepare a magnetically aligned sample. The 
X-ray diffraction data on the aligned sample indicates 
the [010] direction preferably oriented along the magnetic 
field. The degree of alignment is estimated to be 42 ± 
10 % by a method described later. The other crystal 
axis directions may be randomly oriented in the plane 
perpendicular to the magnetic field. 

III. SUSCEPTIBILITY AND NMR 
MEASUREMENTS 

The susceptibility measurement has been made by a 
SQUID magnetometer with an applied magnetic field of 
4 T. Fig. |21 shows the polycrystal susceptibility of the 
epoxy-free sample. Among the literatures of suscepti- 
bility measurement, Ref. [l0| is the only work which has 
reported the high temperature data beyond the suscep- 
tibility maximum. The small discrepancy found between 
(A) and (B) is unknown. We will use below the data of 
Ref. ^3 to make a comparison in the isotropic part of spin 
susceptibility between the experiment and theory. 

We need an anisotropic susceptibility data to get 
an orbital susceptibility which is primarily anisotropic. 
The uniaxially anisotropic susceptibility of the present 
magnetically-aligned sample has been obtained as is 
shown in Fig. 0] Single crystal data by Ref. shows a 
larger anisotropy than the present sample, provably be- 



FIG. 4: Anisotropic susceptibility of CuO. The data set of 
(a) is the magnetically [010] aligned polycrystal susceptibility 
by the present work, the data set of (b) is the single crystal 
susceptibility taken form Ref. llH , and the data set of (c) is 
given by an alignment correction to the data set of (a) by the 
method described in the text. The origin of the offset found 
in the absolute value of (b) is unknown, although we have 
observed a similar offset in a single crystal prepared in our 
Institute fRef. liH). 

cause the present magnetic alignment is only partial. The 
temperature range of the single crystal data by Ref. IT3I 
unfortunately, is limited below 264 K, too narrow to meet 
the following analysis together with NMR data we mea- 
sured. Thus we use two types of susceptibility data in the 
following analysis to get the orbital susceptibility; one is 
the raw data (a) as is measured, another is the corrected 
data (c) given by a method described below. It can be 
seen later that the main conclusion drawn from the fol- 
lowing analysis turns out to be independent of which sus- 
ceptibility data are used. 

The corrected data have been taken as follows. We as- 
sume the present sample is a sum of crystallites oriented 
perfectly and polycrystalites oriented randomly, then we 
can write, 

4=Pxl + l(l-p)(xl+2Xs)> (la) 



Xp =Pxi + ^-p)(xl + ^xi), (lb) 

where p is the degree of alignment, \v anc l Xs are the sus- 
ceptibilities of partially aligned sample and single crystal 
one, respectively, measured in a field parallel (x") and 
perpendicular (x' 1 ) to the [010] direction. Comparing 
the anisotropy of the raw data (a) with that of the sin- 
gle crystal data (b) at 264 K gives p = 0.42 ± 10. The 
degree of alignment estimated in this manner is found 
nearly temperature independent in the paramagnetic re- 
gion, as is expected. We can estimate the fully aligned 
susceptibility in the temperature range of 231 — 400 K by 
substituting the p and raw data to the inverse of Eqs. 
The result is shown in Fig. 0] (c) . 
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FIG. 5: 63 Cu-NMR -1/2 «-> +1/2 spectrum taken at 300 
K and 68.0 MHz in the present sample (not aligned). The 
spectrum shows a powder pattern which undergoes the 2nd 
order quadrupole effect as well as magnetic shift. Each point 
has been taken by accumulating about a thousand of spin echo 
signals. The signal has become hard to be observed below 240 
K, because of the critical slowing dawn of spin fluctuations 
towards the magnetic phase transition at 231 K. The arrows 
indicate the particular resonance positions of which we took 
the frequency dependences to deduce the NMR shifts K. 
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FIG. 6: The Cu-NMR paramagnetic shift with respect to the 
effective gyromagnetic ratio of 63 7/(2vr) = 11.284 MHz/T is 
plotted against temperature. We measured the anisotropic 
shift for the two particular field directions parallel and per- 
pendicular to [010] axis at each temperature. The Tn2 (231 
K) indicates the 2nd Neel temperature. 



Paramagnetic 63 Cu-NMR shift K is obtained by mea- 
suring a frequency dependence of the spectrum (Fig. O, 
and the result is plotted in Fig. In order to de- 
duce the magnetic contribution of K separately from the 
quadrupolar one in the spectrum, we have eliminated the 
quadrupole effect j2jj by taking a frequency dependence 
of the particular resonance positions indicated by arrows 
in FigEl 

The orbital susceptibility is deduced, as is listed in 
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FIG. 7: The Cu-NMR shift is plotted against susceptibility 
with the temperature as implicit parameter. The data set 
of closed symbols (a) are given from the susceptibility raw 
data (a) of Fig. the data set of open symbols (c) are from 
the corrected susceptibility data (c) shown in Fig. |3] respec- 
tively. The solid and dotted lines are linear fits through the 
respective data sets. 



Tab. II I II by a well-known method [21j of analyzing the 
K~x plot diagram (Fig. 0) together with the following re- 
lations between the total shift K and total susceptibility 
X written as, 



and 



K = K spm (T)+K 



X = Xspin(T) + Xorb + Xc 



A s ■ 

K S pi n (T) — — Xspin(T), 



N A 



'Xorbi 



(2a) 



(2b) 



(3a) 



(3b) 



where we use the usual notations as the hyperfine cou- 
pling tensor A sp i n including the Fermi contact and dipole 
interaction, the Avogadro's number N^, the Bohr mag- 
neton \xbi expectation value (r -3 ) within the 3d ra- 
dial wave function and the covalence reduction factor 

o< i 2 < i. 

The orbital susceptibility Xorb tells us the g factor, as 
listed in Table IITT1 by using the following relations as [2I| , 
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TABLE III: The orbital susceptibility Xorb, crystal field split- 
ting A, and g factor are obtained by analyzing the straight 
lines shown in Fig. |7| together with Eqs, 12151 The column (a) 
and (c) are obtained from the data sets of (a) and (c) in Fig.[7| 
respectively. The susceptibility and crystal field splitting are 
in the units of 10 -4 emu/mol and eV, respectively. 
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0.139 
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and 

ff"=2-%£, (5a) 

9 =2 ^— , (5b) 

where we denote the spin-orbit coupling A so of the 3c? 9 
state, and the crystal field splittings A and of Cu 2+ 
ions in a tetragonal crystal field. 

In Table ITTT1 we assume the following values, as com- 
monly expected for cuprates [2^ , that the total ionic core 
susceptibility of Cu 2+ and O 2- is — 0.23xl0 -4 emu/mol 
[23| . the hyperfine radius parameter (r~ 3 ) = 8.25 a.u., 
covalence reduction factor £ 2 = 0.75 [24J and spin-orbit 
coupling A so = -0.095 eV |25|. 



IV. SPIN SUSCEPTIBILITY AND 
SUPEREXCHANGE INTERACTION 

The experimental spin susceptibility Xspin can be ob- 
tained by subtracting the orbital and ionic core suscepti- 
bilities from the observed total susceptibility. For the 
total susceptibility to be subtracted, we use here the 
data of Ref. ^] rather than the present data, in order 
to show the whole behavior including the susceptibility 
maximum. We take the isotropic component of orbital 
susceptibility from Table|m]as xSrt = §(xlrt + 2 Xorb) = 
0.31 x 10 _4 emu/mol and the literature value of Xcore — 
—0.23 x 10 _4 emu/mol [23j for the subtraction. We note 
here that the difference in the orbital susceptibility be- 
tween the column (a) and (c) in Table ITTT1 is about 3 % 
of the observed total susceptibility, being small enough 
to make only a negligible difference in the result of spin 
susceptibility. This is also related with the principle that 
the spin susceptibility of a S = 1/2 system is primarily 
isotropic. 

In F ig. El the dotted curves are theoretical predictions 
for ID |26j and 2D [23 5 = 1/2 Heisenberg antiferromag- 
net. Since the theoretical predictions are given in a scale 
normalized by the coupling constant J for both temper- 
ature and susceptibility axes, we need a particular value 
of J to make a comparison possible between experiment 




kJ/J 

FIG. 8: Comparison between experimental and theoretical 
spin susceptibility. The dotted curves are theoretical predic- 
tions for ID [2(J and 2D S = 1/2 Heisenberg antiferro- 
magnets. We get the experimental curves "ID" and "2D" by 
applying J = 73 ± 3 (ID) and 52 ±3 meV (2D), respectively, 
to the experimental spin susceptibility and temperature. 

and theory. Each theory predicts a relation between J 
and T max given by, 

J = k B T max /d, (6) 

where T max is the temperature showing susceptibility 
maximum, d is the numerical factor given by d =0.64 
and 0.9 for ID and 2D, respectively, and J is defined in 
a spin Hamiltonian expressed by, 

H = -JStS 2 . (7) 

Putting the experimental value T max = 540 ± 20 K [l(| 
into Eq. © we get J = 73 ± 3 and 52 ± 2 meV for ID and 
2D, respectively, and correspondingly the experimental 
curves "ID" and "2D" , respectively, as shown in Fig. [8] 
The g factor used in the normalized susceptibility axis is 
given by substituting the values listed in Tab. lIIII into the 
isotropic component expressed as g = ^(g^+2g ± ) = 2.09. 

We can find from Fig.[H]that the ID theory reproduces 
nicely the experimental ID curve above the Neel tem- 
perature T/V2. The 2D case shows a poorer agreement 
between the 2D theory and the experimental 2D curve 
at all temperatures. This gives a support for the hypoth- 
esis that a quasi ID AFH model is a good approximation 
of CuO. 

We emphasize here that the discrepancy between the 
2D theory and the experimental 2D curve is near the 
isotropic component of orbital susceptibility, which sug- 
gests that estimating orbital susceptibility is crucial to 
make a comparison in the absolute value of spin sus- 
ceptibility. If the orbital susceptibility is taken as an 
adjustable parameter, as is often assumed, it would be 
hard to make a distinction between ID and 2D spin sus- 
ceptibility. 
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TABLE IV: The geometrical parameters of Cu-O-Cu bond 
(bond angle 6 and bond length L) and superexchange cou- 
pling J in the antiferromagnetic cuprates are listed. The J 
is defined by Eq. Qin the text or by the 1st nearest superex- 
change coupling in the case of CuGeOa. The data of J are 
taken from experiments of magnetic susceptibility, neutron 
and Raman scattering measurements. The bond angles and 
bond lengths are obtained by calculations using the litera- 
ture data about the crystal structures; CuGeOs from Ref. l2£l: 
CuO from Ref. [3 YBa 2 Cu 3 6 , La 2 Cu0 4 and Nd 2 Cu 4 from 
Ref.lgl The and L of CuO refer to the bond Bi 22 . 
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FIG. 9: The superexchange coupling J is plotted against the 
bond angle 8 and against the bond length L (the inset) in the 
cuprates. The data are taken form Tabled The solid curve 
shows a cubic fit through the points. The triangles are the 
data taken from neutron scattering, the squares from Raman 
scattering, the circles from susceptibility measurements. 



The earlier result of J = 69 ± 20 meV by neutron mea- 
surement |2( is in better agreement with the present value 
of 73 ± 3 meV (ID), than that of 52±2 meV (2D). This 
yields another support for the quasi ID antiferromag- 
netism of CuO. The present agreement may be sufficient 
to argue that the effect of ID quantum spin fluctuations 
[26] | is responsible for the reduction of the spin moment 
0.65 [Ib compared with 1 /is expected for a Cu 2+ ion. 

The chain direction of [10—1] may be most promising 
from the viewpoint of the bond angle of the bond Bi22- 
The advantage by the bond angle of B122 would, never- 
theless, trade off the disadvantage by the bond length of 
itself. A careful explanation for the roles of bond angle 
and bond length in magnetic coupling J is required. We 
will be focused to this point in the followings. 

Table summarizes the typical examples of J ob- 
tained experimentally for the series of cuprates, together 
with the bond angle and bond length. We plot them in 
Fig- El m which a systematic correlation is found between 
J and 8, but no straightforward correlation between J 
and L. The 8 and L of CuO refer to the bond Bi22- If 
we take the inter-chain bonds B^ft rather than B122, no 
correlation was found in both the plots. 

It can, in principle, be expected that the superex- 
change coupling J would be given by a function of both 
bond angle and bond length. Fig.|§limplies that the bond 
angle is effectively of a prime importance in the magnetic 
coupling of cuprates, and the variable range allowed for 
bond length in the cuprates may be narrower than that 
of bond angle, the latter being capable of a full range 



from 90 degree through 180 degree. 

A bond angle dependence of J in the ferromagnetic 
cuprates has also been found in Ref. m which they 
have reproduced theoretically the bond angle dependence 
of the ferromagnetic J in the vicinity of 9 — 90 degree by 
a close inspection into the contribution from the charge- 
transfer gap A to the superexchange coupling J. This 
tempts us to believe that the bond angle dependence of 
antiferromagnetic J found in the present study is con- 
ceivable. 

We can find from Fig.[§]a reason why CuO is a quasi ID 
antiferromagnet. The 8 dependence of J tells us that all 
the inter-chain bonds Cu-O-Cu, having the bond angles 
ranging from 94.6 to 108.5 degree (see Table |j"T|) . can be 
expected to have the J values smaller than 15 meV. The 
ratios of the inter-chain J values (15 meV or less) with the 
intra-chain J value (73 meV) are 0.2 or less, which may 
be sufficient to make CuO a quasi ID antiferromagnet. 

V. DISCUSSION 

We have assumed during the present study that the 
Heisenberg model is a good approximation for the spin 
system. It could be a possible issue to examine to what 
extent the 2nd nearest neighbor Cu 2+ ion comes into 
magnetic interaction. This may be investigated in a fu- 
ture work, but we can note at the present time that a 
possible 2nd nearest superexchange coupling in the spin 
chain seems negligible, as far as the susceptibility is con- 
cerned. 
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This gives a contrast to the case of the ID antiferro- 
magnet CuGe03 where the spin susceptibility is consid- 
erably reduced, compared to the ID Heisenberg model, 
provably by the Majumdar- Ghosh type of a spin-gap due 
to the significant contribution from the 2nd nearest su- 
perexchange interaction in the chain [^O, The crystal 
structure may be responsible for the difference between 
CuO and CuGeOs. The latter has a edge sharing spin 
chain in which a possible overlap integral between the 
two oxygens locating along the chain may increase the 
2nd nearest superexchange interaction. 

We do not think the bond angle dependence shown in 
Fig. [5] to be an absolute relation, but expect it to be 
an empirical relation valid in some cases. For example, 



considering the mirror and inversion symmetry of the Cu- 
O-Cu bond, the bond angle dependence of J would have a 
zero gradient at 9 — 180 degree, being not consistent with 
the present result in Fig. [5] A broken-symmetry caused 
by the surrounding ions such as alkaline-earth and rare- 
earth metals in Yf^CuaOe and La2Cu04, which are 
located at sites without mirror symmetry on the Cu02 
layer, may be responsible for the singular behavior found 
near 8 = 180 degree. 

One of the authors (TS) expresses his thanks to 
Masashi Hase, Sadamichi Maekawa, Russell E. Walstedt, 
Noriaki Hamada, Takami Tohyama, Hisatoshi Yokoyama 
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